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Abstract. The correlations in classical multi-component ionic mixtures with 
spatial dimension d>2 are studied by using a restricted grand-canonical ensemble 
and the associated hierarchy equations for the correlation functions. Sum rules 
for the first few moments of the two-particle correlation function are derived and 
their dependence on d is established. By varying d continuously near d = 2 it is 
shown how the sum rules for the two-dimensional mixture are related to those for 
mixtures at higher d. 
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1. Introduction 

The statistical equilibrium properties of classical many-particle systems with long- 
range forces have been the subject of an extensive literature (for reviews see [l]-[4]). 
The simplest models with long-range interactions are Coulomb systems consisting of 
point particles with charges of the same sign that move in an inert uniform background 
of opposite sign. For these systems no collapse of particles can occur and stability is 
guaranteed at all densities and temperatures. Both the one-component plasma, also 
known as jellium, and ionic mixtures of particles with different charges and masses 
fall in this class. An important tool in the analysis of the equilibrium behaviour of 
these systems is furnished by the set of correlation functions and the associated Ursell 
functions. The first few moments of the latter satisfy sum rules, which are essential 
for the description of the large-scale fluctuations of local densities. 

In studying one-component plasmas and ionic mixtures it has been found that 
the dimension d of space in which these systems are embedded plays a remarkable 
role. It turns out that several properties of systems with d = 2 and d = 3 (which 
have mainly been considered) are quite similar, whereas occasionally the derivation of 
these properties proceeds along rather different lines. An example is a recent proof of 
a second-moment sum rule for correlations near a guest charge in a two-dimensional 
one-component plasma '5] . Here the use of symmetry properties of the Ursell functions 
leads to a short proof whereas in deriving the analogous sum rule for the three- 
dimensional case a detailed analysis of the statistical ensemble properties has to be 
carried out [7] . Sometimes the analogy between the two- and three-dimensional cases 
gets lost altogether, as seems to be the case for a higher-order sum rule of the two- 
dimensional one-component plasma [8] . For this sixth- moment rule no counterpart at 
d = 3 has been found as yet. 

The purpose of the present paper is to postpone any choice of dimension and to 
derive sum rules that are valid for ionic mixtures in all dimensions d > 2. We shall 
refrain from a discussion of the case d = 1, as periodic oscillations in the density 
lead to complications in that case [2]. We shall concentrate on sum rules for two- 
particle Ursell functions. Our unified treatment enables one to clearly see how the 
simplifications in the derivation of these sum rules for d = 2 come about, and why the 
proof for d > 2 (and hence for d = 3 in particular) is necessarily more complicated. In 
the course of our analysis we shall obtain several new results for a general Coulomb- 
type system with d > 3, which has hardly been discussed in the past [9]- [13]. In 
deriving our results we shall treat d as a continuous variable, as is standard practice 
in the theory of phase transitions [14] and in dimensional regularization of quantum 
field theory [TS] . This method has been used in the context of systems with long-range 
forces as well [16]. 

When describing multi-component ionic mixtures attention has to be paid to a 
suitable choice of the equilibrium ensemble. As in a previous treatment [7 , we shall 
use a restricted grand-canonical ensemble, in which the fluctuating particle numbers 
are constrained by stipulating that the ensuing total charge matches the fixed charge 
of the inert background. 

2. Ionic mixtures in dimension d > 2 

We consider a d-dimensional multi-component ionic mixture of s components, with 
label (T = l,...,s, ina large volume V. The No- particles of species a carry mass 
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and positive charge e„. The system is neutral owing to a uniform background with 
charge density -q^ = ~Ha eaN„/V. 

For arbitrary d the potential (f) depending on the distance r = |r| is proportional 
to l/r'*"^. It is the solution of the d-dimensional inhomogeneous Laplace equation 
A0(r) = — <^(r), with A the d-dimensional Laplace operator and 5{y) the Dirac delta 
function in d dimensions. Here it should be noted that in a space with dimension d the 
Laplace operator acting on an isotropic function is given by [d/dr) r'^~^ [d/dr). 

The explicit form of (f){r) is 

V{d/2 - 1) 1 

'^(^) - 4^rf/2 ^ + (2-1) 

with r(z) the gamma function and with Cd an arbitrary additive constant. We used the 
fact that the surface of a unit sphere in d dimensions equals 2 -k'^/'^ /T{d/2). For d = 3 
the potential has the form (/)(r) — l/(47rr) (at least for Cd = 0), which corresponds to 
the choice of so-called rationalized Lorentz-Heaviside units in electrodynamics. For 
d = 2 the potential that solves the two-dimensional Laplace equation is logarithmic: 

0(r) = log(r) + c (2.2) 

with a constant c that can be used to render the argument of the logarithm 
dimensionless by writing c — log(L)/(27r) with an arbitrary length L. This potential 
can be obtained from (|2.ip by taking the limit d — > 2, if is chosen as : 

r(rf/2 - 1) ^ 

Indeed, in the limit d ^ 2 one finds: 



It should be noted that the shift in energy Cd as given by (|2.3|) becomes infinite, when 
d tends to 2. This does not come as a surprise since the potential (|2.2p grows without 
bound for large r, whereas the potential (|2.ip for d > 2 tends to Cd at large r. If 
desired, one may choose Cd to be given by (|2.3p for all d. However, we shall see that 
for c? > 2 many formulas simplify by choosing Cd — 0, so that the choice p.3p is 
somewhat artificial in that case. For that reason we shall postpone a specific choice 
of Cd and leave it arbitrary as yet. 

The Hamiltonian of the ionic mixture is the sum of the kinetic energy T and the 
potential energy U: 

H = T + U + ^ 6,16.2 '/'(Ir.iai -r^.a. I) 

era ^ cTiai, (720:2 

~qv J2 6- dr (/.(|r,„ - r|) + i ql T dr dr' 0(|r - r'|) (2.5) 

The particle a of species a has position r^Q and momentum Pcra . The prime at the 
summation sign indicates the condition aiai ^ 020.2, so that self interactions among 
the point particles are excluded. As said above, the constant Cd in the potential is left 
arbitrary for the time being. The integrals representing the interactions involving the 
background are taken over the d-dimensional volume V . 

As was shown by Lieb and Narnhofer |17| for the one-component plasma in 
dimension d = 3, the potential energy [/ in if is bounded from below, so that the 
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stability of the system is warranted in that case. Generahzing their argument so as 
to be applicable to a mixture in arbitrary dimension one may prove stability for any 
d > 2, as is shown in appendix A. 



3. Electrostatic sum rules 

The fc-particle equilibrium correlation functions gcn...ak satisfy the BGY hierarchy 



equations [2\ 

d ^ d 

9i\\.aA^u---,rk) = -f3e,, ^jn, . . . , r^) ^ e,^ 0(|ri -r,|) 

rV 

-13 ^ e^,_^, / dvk+i [si'j+.^j.+i (ri, ■ ■ ■ ,rfc+i) - g^^'..^Jri, . . . ,rfc) 

x^<^(|ri-rfe+i|) (3.1) 

with (3 — {kBT)^^ the inverse temperature and n^j ~ {N^y) /V the average particle 
density of species a. The correlation functions can be expanded in terms of Ursell 
functions /iCTi...crfe [13 US]. In particular, the two-particle Ursell function haia2 is 
defined as gaia2 — 1- For large V the Ursell functions are translationally invariant, 
so that they depend on the difference between the positions only. In the following 
we shall assume that the Ursell functions satisfy the standard exponential clustering 
hypothesis, which implies that they tend to zero faster than any power if the separation 
between two positions goes to infinity. 

For k = 2 the hierarchy equation reads in terms of the Ursell fmictions: 



d(t){ri3) _ dh^a}a2 (ri,r2) 



d f 
f^^'^^'O^^'Y^ '^'"36^3 / dvz h^a]as (r2, rs) ^(rig) 



-/^e.,e., e.(ri,r.) ^-/3e.,e., ^ (3.2) 



with rij — Ti — Tj . 

The second term at the right-hand side can be rewritten by expanding the 
potential in terms of Gegenbauer polynomials. For r > r' one has [50]: 

i=i ^ ' 

with d> 2. Here 6 is the angle between r and r'. For r < r' a similar expansion holds, 
with r and r' interchanged. By expanding the potential in this way and using the 
orthogonality relation of the Gegenbauer polynomials one may establish the identity 

_d_ 



(3.4) 



r23<I'12 



Employing this equality in (13. 2p and making use of the exponential clustering 
properties of the Ursell functions one proves the perfect-screening condition [9l [TOl [TTl 
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CT2 



Similarly, by using the Gegenbauer expansion and the exponential clustering 
property one derives from the hierarchy equations for fc = 3 the perfect screening 
rules: 

^n^3 / drg hi%^„^{ri,r2,r3) = -{e^, + e^J h'^^^^{ri,r2) (3.6) 

-~^^^e..ri2h^^lM^V2) (^=1,2,...) (3.7) 

with 9 the angle between Vi2 and ri3, and with {x)n — x{x + 1) • • • (x + n — 1) the 
Pochhammer symbol. 

For d — ?> the Gegenbauer polynomials in (13. 7|) reduce to Legendre polynomials, so 
that we recover one of the well-known perfect-screening rules for a three-dimensional 
ionic mixture [71 [TTJ [53] . To derive the analogous identity for d = 2 we use for £ > 1 

V,va^Cf{x) ^jT,{x) (3.8) 

with Ti{x) the Chebyshev polynomials of the first kind. With the help of this relation 
one finds from p.3p in the limit d 2 the standard expansion of the logarithmic 
potential [5D]. The perfect-screening rule (|3.7p becomes upon taking the Hmit d ^ 2: 



'^'^3 J dr^ (ri , r2, ra) rfg (cos ( 



0-3 

= -e.,rf2/iS2),,(ri,r2) (€=1,2,...) (3.9) 
which for the one-component case corroborates a previous result [H] . 

From the above results a consistency relation can be obtained. On one hand, we 
can prove from (j3.2[) with (|3.4p and (|3.5p . upon multiplying by ri2 and integrating 
over r2- 

d(l){ri3) 



/ dr2 drs h'-^^^^^^ (ri , ra , rg) 

0'3 

= -5/3 e^, ^n„3 e^3 ^ drg /i(,2)^3 (rj, rg) 



23 



/ dr2 4%^(ri,r2) 0(ri2) 
+ [d - (d - 2) Cd /? e,, e^J y dr2 /^^'^)<,^(ri,r2) (3.10) 



On the other hand, from (|3.7p for € = 1 one gets after multiplication by r^2 ^^^d 
integration over r2'. 

^n^3 e^3 / dr2 drg /i(,%^^3 (ri, r2, rg) rig ■ = 

<^3 

= (d-2)e,,y" dv2h'^;^l^{v^,V2)4>{ri2)-{d-2)cde„, j dr2 ^^^.^ (ri, rs) (3.11) 
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Comparison of p.lOp and p. lip yields an identity, which by means of p.5p gets the 
simple form: 

/2 d 
cfr2 ^if^2(ri,r2) r^2 — jr (3-12) 

For d = 3 this identity reduces to the well-known sum rule that was first obtained by 
Stillinger and Lovett [55] and discussed subsequently extensively [71 [121 [HI [HI [571 HB] ■ 
For the one-component case with c? > 3 its form has been found before [T5] . 

The sum rule p.l2p is independent of Cd, as it should be, since the correlation 
functions cannot depend on the choice of an additive constant in the potential. For 
d>2 the intermediate steps in deriving p.l2p simplify for the choice Cd — Q, but that 
is not essential for the proof. To treat the limit d — > 2 one has to choose the specific 
value (|2.3p for Cd, so that (^(r) stays finite. With that particular choice the proof of 
(|3.12p remains valid in the limit d — > 2. The form of p.l2p for 0? ^ 2 is consistent with 
that found previously by taking d = 2 from the start ^24, 28J. The above derivation 
shows how the general form of the Stillinger-Lovett relation for an ionic mixture reads 
for arbitrary d > 2. 



4. Equilibrium ensemble and thermodynamics 

To prepare the ground for the derivation of additional sum rules for the pair correlation 
functions of the ionic mixture we need to specify the equilibrium ensemble for the 
system. A convenient choice, which has been discussed before [7], is the restricted 
grand-canonical ensemble. It is a grand-canonical ensemble with particle numbers 
satisfying the constraint Co- = V . Its partition function Z depends on the 
volume V, the inverse temperature /3, the background charge density g^, and s — 1 
chemical potentials /io-((T 7^ 1). In the limit of an infinite system the partition function 
leads to a thermodynamic function p that is defined by writing: 

lim i logZ(/3,{/3/i4,g„,y)=/3p(/3,{/3A.},g.) (4.1) 

V — ^oo V 

The energy density Uy and the particle densities for a ^ \ follow by taking 
derivatives: 

n.--^ ' "- = 9^ (- = 2,...,.) (4.2) 

In writing a partial derivative with respect to one of the variables /?, 1/3/1^}, (/t,, the 
other variables that are meant to remain constant are suppressed. The pressure follows 
from p through the relation 

p = p-qv (4.3) 

oqv 

as is proved in appendix A. 

For d > 2 a, scaling argument can be used to relate the partial derivatives of p. 
In fact, the potential energy satisfies the identity 

U{r''\...,v''%V)^X^-^U{\v^\...,\r''%X''V) + \ {X^-^^l) CdY^N.el 

(4.4) 
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for arbitrary positive A. This property implies a specific scaling behaviour of the 
partition function Z and the thermodynamic function p. As a consequence, the 
pressure and the energy density of the ionic mixture are related as 

Uy - ^{d-A) - + Cd2_^n„ (4.5) 

a 

with n — rio- the total particle density. 

The partition function Z, and hence p, depends on the additive constant Cd via the 
Hamiltonian. However, the combination (|4.3p . which gives the pressure p, is invariant 
when Cd is modified. On the other hand, the energy density Uy as given by (|4.2p does 
depend on Cd- Its dependence is such that m„ + i rio- is invariant, so that 

(|4.5p can be satisfied. The specific amount by which the Hamiltonian is shifted when 
a different choice for Cd is made depends on the particle numbers N^, as (j2.5p shows. 
Hence, the chemical potentials fl„ (with cr ^ 1) change as well when a different value 
for Cd is chosen. However, the combination fia + \ Cd e^ {e^ — e\) is found to be 
invariant. Of course, the partial densities n„ do not depend on Cd- 

It should be noted that both the pressure and the energy density can be written 
as a sum of a kinetic and a potential part: 

According to (|4.5p the potential parts of the pressure and the energy density are 
related as p^"* — \{d — 2)/d] (u^"' + | n^. e^). In appendix A it is shown how 

several auxiliary relations can be derived from (|4.5p . 

For dimension > 2 one may take Cd — 0, so that (|4.5p gets a simpler form |12j . 
In contrast, for d ^ 2 one should choose Cd according to (|2.3p . With that choice the 
energy density Uy stays finite for d ^ 2. Hence, it drops out from (j4.5p in the limit. 
As a consequence, we are left with the equation of state for the two-dimensional ionic 
mixture: 

"-i-sE"-' (4-7) 

(7 

which can also be obtained directly by applying a scaling argument to a system with 
a logarithmic potential [29l [30] . The present derivation shows how the second term 
at the right-hand side comes about as a consequence of the shift Cd in the potential. 
Incidentally, we remark that it is essential to choose the right value for Cd before taking 
the limit d — > 2. For instance, choosing Cd = in (j4.5p and taking the limit naively, 
without realizing that Uy diverges in that case, would have resulted in an incorrect 
equation of state. 

In closing this section we remark that alternatively one may choose to describe 
the equilibrium ionic mixture by means of a full grand-canonical ensemble with a 
background with fixed density fT2 l [21 1 [3T 1 [32] . 

5. Thermodynamic sum rules for pair correlation functions: zeroth- and 
second-moment rules 

In the restricted grand-canonical ensemble the derivative of the partial density 
with respect to the chemical potential combination Pjla-2 is given by 
Dn^. 1 . 1 . . . . 
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with the operator D/D(3jic! defined as 



a ^1 



The right-hand side of (15. ip can be expressed as an integral over the pair correlation 
function. As a result one finds: 

Ua, J dr2 ft.^^^^^^ (ri , Ta ) = - S„^„^ (5.3) 

Upon summation over (T2 , with the weights Co-j , one recovers the perfect-screening rule 
(|3.5p . Taking an unweighted sum over a2 and using (|A.8|) to eliminate the derivative 
of the particle density n we find the equality: 



id (d - 4) / h^aL (i-i, 1-2) = -(d - 2) /3 ^ + d q„ 



9g« 



-id (d - 2) n., + i (d - 2) /? (^X] 71., e^^ j (5.4) 

This zeroth-moment sum rule is independent of the perfect-screening sum rule. Like 
that rule it is valid for each species cti separately. If an unweighted sum over ai is 
carried out, one arrives at a less strong sum rule of the form: 

i d (d - 4) ^ "^1 J dv2 h'^;^l^{Yi,Y2) = -{d~2)p^ + dq,^ 

-id(d-2)n+i (d-2)c,/3 ^e^^ (5.5) 

For any d > 2 one may choose = in (|5.4p and (|5.5p . For d = 4 the integrals 
in (|5.4p and (|5.5p drop out; the resulting equalities are trivial consequences of the 
relation (jA.Sp . The case d = 2 deserves special attention, and will be discussed at the 
end of this section. 

The derivative of the partial density ricr with respect to the inverse temperature 
(3 reads 

8n 1 1 

-^ = ^-{N,H) + -{N„){H) (5.6) 

Like the derivative with respect to the chemical potentials discussed above, it can be 
written in terms of integrals over Ursell functions, as shown in appendix B: 

{d-2) (3 ^ -I i3q^ n^^ ^ n., e., J dra /i(,^^)<,Jri, ra) rJa 

~\d{d~ 4) ^ I dv2 h^i^}^^ {^1,^2) - \ d {d ~ 2) 

+ \ (d-2) CdPn„, 



(5.7) 



Yn„2 / dY2 h^^}„^{ri,V2) + el^ 

An essential role in the proof of this identity is played by the symmetry properties of 
the Ursell functions, as is discussed in appendix C. Employing (|5.3p and (|5.4p for two 
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of the integrals at the right-hand side, we find that many terms cancel. In this way 
we obtain the second-moment sum rule: 



..|dr,;^(f^,(ri,r,)r?,^-^^ (5.8) 



Summing over ui with the weights e^i one recovers the Stillinger-Lovett rule ()3.12|) . 
Taking the sum with equal weights we get the second-moment identity 



^ ^ ^CJl ^CJ2 ^(T^ 

a"l,(T2 



which is independent of the Stillinger-Lovett rule. 

The above sum rules have been derived for all d > 2. To obtain the corresponding 
rules for the case d = 2 we choose Cd according to (|2.3p and take the limit d ^ 2. The 
zeroth-order sum rule (|5.3p retains the same form, whereas the sum rules (|5.4p and 
(|5.5p become 



^2 

and 



(71 ,CT2 



As in the previous section, incorrect results would have been obtained from (j5.4p and 
(|5.5p when the choice = had been made before evaluating the limit d — > 2. In 
contrast, the sum rules (|5.8p and (|5.9p are independent of the choice of c^, so that the 
proof of their validity for d = 2 is straightforward. It may be noted that in deriving 
the limiting form of the auxiliary relation (j5.7p it is important once again to choose 
Cd correctly before taking the limit. 

The above derivation of ()5.8p for general d shows how one can combine perfect 
screening, symmetry and thermodynamics with the statistical relation ()5.6p to 
establish a second-moment sum rule. For the special case d = 2 the last mentioned 
ingredient is not necessary, as is shown in detail in appendix C. This particular feature 
of the second- moment sum rule (jS.Sp for d ~ 2 has been discovered recently [6] . 

6. Thermodynamic sum rules for pair correlation functions: 
fourth-moment rule 

To derive an equality for the fourth moment of the two-particle Ursell function we 
start from an expression for its derivative with respect to the inverse temperature: 

d r 

- 2) /3 7^ na, hfl^{vi,V2) 



7-2 
' 9 



d(3 

= -I [3 ^n^^ e^3 / drg h^^}^^^^{Yi,V2,Vj,) 

/d 
/iif^,^3(ri,r2,r3) - n^, n^, ri2 • — /i^^^^^Jri, ra) 

-i/3g„n^i /i(,^^)^^(ri,r2) - d (d - 2) /i^^^^^^ (ri, ra) 
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^'^3 ^^3 / *3 ^tf^2<T3 (l-l, 1-2, rs) 



(6.1) 



The proof of this identity is sketched in appendix B. Multiplying both sides with 
Ccri Gcts 4>iri2), integrating over r2 and summing over ai and cr2, we get an expression 
for the derivative of the potential-energy density (jB.ip : 



id(d-4) ^ 



^CTi ^(72 ^(73 ^CTi 6(72 ^CTS 



^CTi ^(72 ^0*3 6fT 



fTi ^(72 ^CTi 6(T2 



dr2 ri 



J dv2 dry, /i(,f^^^^(ri,r2,r3) r^g (/)(ri2) 
, e^, y dr2 dr3 /i^.^^^^^^^ (ri, rj, rg) (/)(ri2) 
0(r-i2) 



9^4S.(ri,r2) 



4 /3 g„ ^ n^, n^, e^^ e„, J dr2 4f^^(ri,r2) 0(ri2) - 2 d (d - 2) <°* 



-i (d - 2) c<j /? ^ ncri n^rj Bo-i Bo 



+ « + y '^'^2 /iL',U(ri,r2) </.(ri2) (6.2) 

The first term at the right-hand side can be expressed in moments of the two- 
particle Ursell functions by using the relation (jC.4p , which follows from the symmetry 
properties of the three-particle Ursell function. Likewise, the second and the sixth 
terms can be rewritten by means of the symmetry relation (jC.2p . In the third term 
we can carry out a partial integration and use the identity 

d(j){ri2) 



ri2 • 



dri 



id -2) cj){ri2) + {d-2) Cd 



(6.3) 



As a result we arrive at a relation involving the zeroth, the second and the fourth 
moments of the two-particle Ursell function: 

+2d{d- 6) {d + 2)l3q^Y^ / dr2 h'^^„^{ri,r2) rl2 

CTl,(T2 

-4 d^ {d - 4) (d + 2) ^ J dv2 /i^'U (ri, r2) = 



{d-2f {d + 2) 0- 



+2 {d-2) {d + 2) Cdf3 



2 '-'"■v 



d(3 



16 {d - 1) (d - 2) {d + 2) (3 ul°' 



(l"l,I"2) 
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-d {d-6) 



7 , ^cri e^^ 



^2{d-2f (d + 2) cl 



2 2 



(6.4) 



The zeroth and second moments at both sides of this relation can be replaced by 
the thermodynamic expressions given in (j5.3p and (jS.Sp . The ensuing derivatives 
with respect to the chemical potentials may be eliminated with the help of (jA.8|) . 
Furthermore, at the right-hand side the full energy density can be introduced with 
the help of (14. 6p . These manipulations lead to an expression for the fourth moment 
of the Ursell function in terms of thermodynamic derivatives only: 



d ql 



E 



dr- 



4!U(ri,r2) 



= -8 (d ■ 
+4 d (d - 



2)2 (d + 2) 



16 (d- 1) (d-2) (d + 2)/3■ 



2) (d - 4) (d + 2) /3 — + 4 d^ (d - 4) (d + 2) (J, — 
op 



-4 (d-2) (d + 2) 



(d-2)/? 




(6.5) 



d , d 
op dqy 

Upon using (jA.9|) we find that the right-hand side is proportional to the derivative of 
the pressure p (in the form of l|4.5p ) with respect to g^,. In this way we have found the 
rather elegant fourth-moment sum rule 

n-cTj riry^ err, errn f dY2 /il^' ''"^ v.^^ ^4 _ 8d(d + 2) dp 

a\ ,£72 



^172 ^<y\ ^<T2 



„i^,(ri,r2) rj^ 



(3qy dqy 



(6.6) 



It is a generalization to arbitrary s and d of the well-known compressibility rule that 
has been established for the one-component plasma in two [H] and three [1111331 [MUSS] 
dimensions and for the three-dimensional ionic mixture [3 HI] . Whereas the second 
moments, as given by the sum rules (I3.12p . (|5.8p and (|5.9p . are linear in d, the fourth 
moment turns out to be quadratic in d. 

For all d > 2 we may put Cd — 0, as before, so that the relations (|6.ip - (|6.5p 
become somewhat simpler. To discuss the case d — 2 we must choose Cd as in (|2.3p . 
Upon taking the limit d — > 2 the terms in (|6.ip ~ (|6.5p containing Cd remain finite, 
so that they cannot be omitted. However, the final result (|6.6p does not depend on 
Cd explicitly, so that it remains valid as such in the limit d — > 2. Hence, we have 
established the fourth- moment rule (|6.6p for all d > 2. 

As a final remark we point out that a shorter proof of the fourth-moment rule for 
the special case d = 2 can be found from a particular symmetry relation connecting 

f 21 

second and fourth moments, as discussed in appendix C. The derivative dhaia2/d(3 is 
not needed in that line of reasoning. 



7. Concluding remarks 

By making a systematic use of the properties of the restricted grand-canonical 
ensemble and the hierarchy equations for the correlation functions we have been able 
to derive the sum rules that govern the first few moments of the two-particle Ursell 
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functions for a multi-component ionic mixture with an arbitrary spatial dimension 
d > 2. The dependence on d of the various moments has been determined in detail. 
While most discussions in the literature had to treat two-dimensional mixtures with 
a logarithmic potential as a separate case, we have shown that a unified description 
of mixtures for all d > 2 is indeed possible by making a careful choice of additive 
constants in the potential. 

Our main results for the moments of the two-particle Ursell function are presented 
in ([331), (|3?T2| . dSSll-dES]), jElll-jSJl) and ([SH). The ensuing results for the moments 

(2) . (2) 

of the two-particle correlation function follow by replacing /io-/ct2 with gaia2 ~ 1- 
Whereas the zeroth-moment perfect-screening rules p.Sp and the second-moment rule 
(j3.12p could be derived without invoking thermodynamical properties, the proof of 
the other sum rules had to be based on statistical ensemble theory. Accordingly, 
the ensuing rules in sections 5 and 6 depend on thermodynamical derivatives with 
respect to the basic variables describing ionic mixtures in a restricted grand-canonical 
ensemble, viz. /3, {Pjla} (for cr = 2, . . . , s) and q„. 

If one wishes, one may express the sum rules in terms of derivatives with respect 
to a different set of independent variables involving - apart from fi - the chemical 
potentials {/io-} with cr = 1, . . . , s, in a way described previously "T. For completeness 
we give the sum rules (15. 3p . (|5.8p and (|6.6p in terms of derivatives with respect to 
these alternative variables: 

^Kfe 4 It: fe) 



^ ^ fT'O'i ^£72 ^(Ti ^(72 



/rfr.e,(r„r.).^, = -^^|±^ (7.3) 



with the abbreviation S = '^^a Gadq^ / djjia- As before, in writing the partial derivatives 
at the right-hand sides the independent variables that are kept constant are suppressed. 

The sum rules discussed in this paper are essential in understanding the 
equilibrium fluctuations in an ionic mixture. In particular, the fluctuations in the 
partial densities, the pressure and the energy density are governed by these rules, as 
has been shown in [36] for the three-dimensional case. The fluctuation formulas in 
turn are necessary in order to determine specific dynamical properties of the ionic 
mixture, such as the time evolution of the collective modes [37[[38l. 
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Appendix A. Stability, thermodynamic pressure and some auxiliary 
relations 

In this appendix we shall first discuss the stability of ionic mixtures in arbitrary 
dimension. Furthermore, we shall establish the relation between the thermodynamic 
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pressure and the partition function in the restricted grand-canonical ensemble. Finally, 
a few thermodynamic auxiliary relations will be derived. 

By generalizing the argument given in [17j so as to be applicable to a mixture in 
arbitrary dimension d > 2 one finds the bound: 

which for the one-component case (and Cd — 0) agrees with the bound presented by 
Sari et al [321 HO] ■ Taking moreover d — 3 one recovers the result in [T7] . For d — > 2 
and Q = the bound in (|A.1[) goes to — oo , so that it becomes useless. However, 
upon choosing q as in (|2.3p the inequality (|A.1[) becomes in the limit d ^ 2: 



U > 



For the one-component case (and the choice c = 0) this inequality has been derived 
previously [39j . It should be remarked that different bounds have been obtained in 
the past [1T]-[11]. For our present discussion these are not relevant, since we only wish 
to confirm here that the multi-component ionic mixture is stable for arbitrary d > 2. 

Furthermore, we want to derive the relation (|4.3p between the the pressure p and 
the thermodynamic function p, which follows from the partition function according to 
(j4.ip . Generalizing the definition of the (thermal) pressure in a one-component plasma 
by Choquard et al [45] to an ionic mixture, we write it as the derivative of the free 
energy F with respect to the volume V at constant temperature T, constant (average) 
particle numbers ricrV (for a = 2, . . . , s) and constant total background charge QvV: 

V "1/ /T,{«„y},?„v 

Taking account of the implicit dependence on V we get 

p=-/.+ ^ fl^) +qJ^) (A.4) 

with fv{T, {ug-}, Qv) = F/V the free energy density. The construction of the restricted 
grand-canonical ensemble implies the relations 0: 



dfv = -s,„ dT + fl^ dn„ + fig dq^ (A. 5) 

<t(#1 



P = -fv + ^ /i<T ric, (A. 6) 

with Sy the entropy density and jig — —dp/dq^. Hence, (|A.4[) can be written as: 

P = - fv + ^ fia n-a + fig qv (A. 7) 

Comparison of (jA.6[) and (jA.7[) yields the relation between p and the pressure p that 
we wished to prove. 

Finally, in the main text we need several equalities involving partial derivatives 
of thermodynamic quantities. Upon differentiating the relation (|4.5p with respect to 
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/3/icr, at constant /3 and q^, we get: 

id-2) P = dq^ ~dn„~\d{d-A) 



with the operator D/Dfijlcr defined in (|5.2p . Likewise, differentiation of (|4.5p with 
respect to fi yields: 



(A.9) 



For d > 2 we may choose Cd — 0, so that the last terms at the right-hand sides 
of (jA.Sp and (jA.9[) drop out. For d —> 2 we choose Cd as in (|2.3p . When the limit is 
taken, the left-hand sides of (jA.8|) and (jA.9[) disappear, while the last terms at the 
right-hand sides yield a finite contribution. As a result we get for d = 2: 



(Iv - n„ + -f^^ \ n - ^^n„, el, \ = (A.IO) 



dn„ Dip 
dqy " Df3fl„ y 8n ^ 

and 

The auxiliary relations (jA.Sp and (|A.9p have been used in the main text. 

Appendix B. Derivatives of densities and Ursell functions with respect to 
the inverse temperature 

In deriving the second- and fourth-moment sum rules we need expressions for the 
derivatives of the partial densities and the two-particle Ursell functions with respect 
to /3. The derivative of n^r with respect to (3 follows by evaluating its formal 
expression: dn„/dl3 = -{N„H) /V+{N„){H) /V . The average (H) of the Hamiltonian 
is proportional to the internal energy Uy , which is the sum of a kinetic and a potential 
part of the form (|4.6p . The latter can be written as an integral over the two-particle 
Ursell function: 

i E "'^1 "'^2 Cai 6,^2 J dr2 h'^^^^^{ri,r2) (j){ri2) (B.l) 

(Tl ,0-2 

The average {NcH) can likewise be expressed in terms of integrals over Ursell 
functions. As a result we find 



d(3 



- = E "'^2 "<T3 6^2 j dv2 drs /i^^^^^^^^ (ri, r2, ra) 0(r23) 



^2) <l>{ri2) 

-\dn„^^ n<,2 y dv2 hS^l^^ (ri, rz) - i d (B.2) 
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Employing the symmetry relation (IC.2p to eliminate the integral with the three- 
particle Ursell function, we arrive at (15. 7p . 

Furthermore, we need an expression for the derivative of the two-particle Ursell 
function with respect to (3. In the restricted grand-canonical ensemble one has quite 
generally d{f)/df3 = — (/ H) + (/) (H) for an arbitrary phase function /. Taking 
/ = ~ ^Tiai) ^(r2 - i"^2a2) and using one derives 

4 



''fT3 "'(74 '--CTs '■^0'4 



J drs dr4 hi%^^,^^^{ri,r2,r3,r4) 0(r34) 
-\dn„, n„,_^na^ j drg /i^.^^^^^ (ri, r2, rg) 

-/3 ri<^2 ^n^3 e,^., ^ drg /i(f^)^^^3(ri,r2,r3) [e^, (^(ris) -f e^^ '^'(''23)] 



-d /i(,^^)^^(ri,r2) ~ P r 

~ 2 P ^(Ti ^CJ2 ^ ^ ^0"3 ^(T4 6a"3 ^(T4 



''0-2 ^CTi ^(T2 "'C7icr2 



/iS^) (ri,r2)0(ri2) 



/i(f,3(ri,r3)/i(2),,(r2,r4) 



(ri,r4)/i(2U(r2,r3)l 0(r34) 



-/? ria, n„.^ ^ i 



dr3 



e<Ti /ii-^,l,(r2,r3) 0(ri3) -t- e^^ ^S^l, ''s) '/'(?-23) 



-/3 rio-i '^cr2 e(Ti 60-2 4>{ri2) 



(B.3) 



For large separation of the position arguments the left-hand side vanishes faster than 

any inverse power of ri2. At the right-hand side, the first three integrals and the two 

(2) 

terms proportional to /ict/(T2 share this feature. However, the property of being of short 
range is not obviously true for the last two integral terms, while it is certainly false for 
the final term, which is proportional to (j){ri2) and hence of long range. Nevertheless, 
by employing (|3.4p and (j3.5p one may rewrite the sum of these terms in a form that 
shows their short-range character as a function of ri2 explicitly: 



d-2 



0'3,0'4 

X 



''(73 ' ^(^S '^(74 



J drs /i^%3(ri,r3) ri3 



d(l){r2?,) 
dr2 



+(3 n^ri ?lcr2 E 



H~/5 ^CTi ^(72 ^(72 



/ /lit, 

-'r24>r23 
^(T3 ^(T4 ^^0-3 ^(74 

X / c?r4 /ii2),^(ri,r4) 

[0(ri2) - Cd] ^ n<T3 / drz /itf^3 (ri , r3) 



(B.4) 



Substituting these terms and using moreover the symmetry relation (|C.6p in the first 
term at the right-hand side of (jB.Sp . we arrive at the somewhat simpler expression 
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(|6.ip given in the main text. It should be noted that at the right-hand side of (|6.ip 
the potential does not occur explicitly any more. 

Appendix C. Symmetry relations 

The Ursell functions are symmetric under a permutation of both their position 
arguments r^ and their component labels <Ji. From that symmetry one proves 

(d - 2) ^ n^, e„^ e^,, j dv2 drg h^^}^^^^{vi,Y2,v^) [(j){r2?,) - c<j] = 

dT2 drs /i(,f^^^3(ri,r2,r3) ri2 • ^^^^^ (C.l) 



At the right-hand side we use p.lO[) . Employing moreover the perfect-screening 
relations (|3.5p and (|3.6p we get the identity: 

T^(J2 '^0"3 6o-2 60-3 / dra dra /i^-f^^^^ (ri, ra, rg) ^(raa) 
+2 {d-2) fj e^, n^, e^^ / dra /i^-^^V^ (^i, ^a) ^(ria) = 
= /5 9^ X! '^'^2 e<T2 y" dr2 /i^!^^, (r 1 , r2 ) - 2 d ^ n^, ^ dr2 ft^^^^^ (r 1 , r2 ) 



J2 "^2 e', y dr2 h^^^,^ (n, rz) + e; 



(C.2) 



0-2 

which is used in section 6 and appendix B. 

A second identity is obtained by starting from an equality that is analogous to 
(jC.ip and follows likewise from the symmetry of the three-particle Ursell function: 

71(j2 ^^(73 ^(T2 ^0"3 

/ dr2 drg /i(,3^)^^^3(ri,r2,r3) [(t){r23) - Cd] 

X [r?2 + ri2 • r23 + 2 r^^ (ri2 • r23)^] = 

dr2 dra /i*,^^)^^^^ (ri, r2, rg) r^2 i"i2 • g^^'^ (C.3) 

cr2,cr3 

The right-hand side can be expressed in terms of two-particle Ursell functions by 
employing the hierarchy equation (j3.2p . At the left-hand side we may invoke the 
perfect-screening rule (|3.7p for € = 1,2, when the sum over ai with weights n^-^ is 
carried out as well. In this way we arrive at the identity: 

{d-2) P ^ n^^ e^^ j dv2 drs /i^.^^^^^ (ri, r2, rg) r^^ 0(^12) 

+ {d -2)(3J2 "'^1 "'^2 e,,, y dr2 /i^^^^^Jri, rz) ^(^2) = 
= 2(^/^2) H ""^1 '^'^2 e^, y dr2 h'^^l^{ri,r2) r\^ 
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-2d ^ e^i J dr2 /i^^^^^^ (ri, ri) rj^ 

-{d - 2) Cd /? ^ el^ e^^ J dr2 h'^^^^^{ri,r2) rl^ (C.4) 



which is needed in section 6 of the main text. 

Finally, we want to establish an equality for the four-particle Ursell function. It 
follows by starting from an equality for ft,'^'*' of a similar form as (jC.ip : 

(d - 2) n^, e<,3 e^, j dr^ dr4 /i^.^^^^^^^^Jri, ra, rg, [0(r34) - Cd] = 

<''3,0'4 

= 2 ^ ri<T3 ria^ e<,3 e^, y drg dr4 /i^-f^^^^^^ (ri, ra, rg, r4) rig • (C.5) 

Upon using the hierarchy equation (|3.ip for fc = 3, the expansion p.Sp . the identity 
p.4p and the perfect-screening rules p.5p ~ p.6p we get, by taking steps analogous to 
those of appendix B of ref. [7]: 

'^(Ts ^(74 6(73 G(T4 

/ drg dr4 /i^.'*^)^^^,,^^ (ri, rs, rg, r4) (/)(rg4) = 

<^3 ,"'4 

f^(T3 60-3 / dr^ ft.^^-'g.^^^ (ri , r2, rg) [eo-j 0(r-ig) + ^(''23)] 
/ f^i"3 ^i^il2CT3(ri,r2,r3) ?-23 

0-3 

c?r3 /iif^,^3(ri,r2,rg) -|-ri2 • — h^^^^^{ri,r2) 

"3 

-P "'^3 "-^-i ^'^3 e<T4 y drg /i(2)^^(ri,rg) rig • / rfr4 /ii^,'^Jr2, r4) 

(r2,i"3) ['/'(^is) - q] 

£'"3, "■4 

X / drih^^^„^{ri,ri) 

Jri4>ri3 

+ {d - 2) /3 [^(''12) - Cd] Y ""^a ^'^3 / '^I'g /i(.2)^3 (ri, rg) 

CT3 •^I-13>''12 

-{d-2) /3e^, e^, h';^^l^{ri, r2) (pir^) + ^ l3 e^, h';^^^^{ri,r2) rf^ 
-\{d-2)cdp ^"-3^3 / drg/it3),^,3(ri,r2,rg) + (e2^+e^J /iS'U(ri,r2) 

. £^3 

(C.6) 

This rather complicated identity has been used in appendix B. Inspection of the terms 
at the right-hand side shows that for large ri2 each of these vanishes faster than any 
inverse power of ri2, as it should be in view of the short-range character of the four- 
point Ursell function at the left-hand side. 

For d > 2 the above identities may be simplified by putting Cd equal to 0. That 
choice is not allowed when one is interested in the limit d ^ 2. In that case one 
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takes Cd according to (|2.3p . In the limit d ^ 2 the left-hand side of the identity (|C.2p 
vanishes, so that we get an identity that connects the zeroth and second moments of 
the two-particle Ursell function: 

P qvJ2'^''2 e^2 J dv2 /iS2^)^^(ri,r2) rl^ = 4^n<,, j dra /i^^^)^^ (ri, ra) 



2tt 



^n„, e% / dr2 /i(.%,(ri,r2) + e^^ 

. CT2 



(C.7) 



It should be noted that the last two terms would have been missed when in (|C.2p 
the limit d ^ 2 had been taken naively after putting Cd — 0. The identity (|C.7[) . 
which is valid for the special case d = 2 only, has been obtained recently [6J. Upon 
substituting (|5.3p in the right-hand side and using (|A.10[) we recover (|5.8p for d ~ 2. 
In fact, this shows that for d — 2 the second-moment sum rule (I5.8P can be derived 
from perfect screening, symmetry and thermodynamics alone, without having recourse 
to the rather complicated expression for the derivative dna/d(3 of the partial density 
with respect to the inverse temperature. The latter expression is essential in deriving 
the second-moment sum rule for arbitrary d > 2. 

Similarly, for d ^ 2 the symmetry relation (jC.4[) reduces to an identity connecting 
the second and fourth moments of the two-particle Ursell function: 

/3 g„ ^ n„2 e^i j dr2 /i^f^^ (ri, rz) r^^ = 

o-i,cr2 

= 16 ^ n^, J dr2 /i^^^^^Jri, ra) 

0-i,(T2 

2(3 
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el^ e^^ j dv2 hf}„^{vi,V2) rj^ (C.8) 

o"l ,a"2 

The last term is missed when one puts = in (|C.4p before taking the limit d ^ 2. 
Substituting (|5.8p and using the equation of state (|4.7p we are led to (|6.6p for d — 2. 
Hence, in a similar way as discussed above for the second-moment sum rule, the 
derivation of the fourth- moment sum rule can be simplified for the special case d ^ 2. 
For that case it is enough to make use of the perfect-screening and second-moment 
rules, symmetry properties and thermodynamical relations in the proof, whereas for 

(2) 

general d the derivative 9/io-i 0-2 / dP of the two-particle Ursell function with respect to 
the inverse temperature needs to be determined. Incidentally, we remark that for the 
one-component case the identity (|C.8P has been obtained before [24] . 

We are left with (|C.6p in the limit d ^ 2. The resulting identity is rather 
complicated and is not needed in the main text, so that we refrain from writing it 
down. 
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